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ABSTRACT 

We investigate the tidal interactions of a red giant with a main sequence in 
the dense stellar core of globular clusters by Smoothed Particle Hydrodynamics 
method. Two models of 0.8 M red giant with the surface radii 20 and 85R & 
are used with 0.6 or O.8M main sequence star treated as a point mass. We 
demonstrate that even for the wide encounters that two stars fly apart, the 
angular momentum of orbital motion can be deposited into the red giant envelope 
to such an extent as to trigger rotational mixing and to explain the fast rotation 
observed for the horizontal branch stars, and also that sufficient mass can be 
accreted on the main sequence stars to disguise their surface convective zone 
with the matter from the red giant envelope. On the basis of the present results, 
we discuss the parameter dependence of these transfer characteristics with non- 
linear effects taken into account, and derive fitting formulae to give the amounts of 
energy and angular momentum deposited into the red giant and of mass accreted 
onto the perturber as functions of stellar parameters and the impact parameter 
of encounter. These formulae are applicable to the encounters not only of the 
red giants but also of the main sequence stars, and useful in the study of the 
evolution of stellar systems with the star-star interactions taken into account. 

Subject headings: globular clusters: general — stellar encounters: tidal captures 

1. INTRODUCTION 

In the core of globular clusters, it is thought that star-star interaction s play an importan t 



np 

m 



role because of very high stellar density and of relatively low velocity (IHills &: Daylll976l ). 
There is growing evidence for the modifications of stellar properties and populations under 
the influence of close encounters and collisions. For example, the overabundance of low-mass 
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X-ray binaries and millisecond pulsars is regarded as consequent upon tidal captures of an 



Fabian et al. 


1975; 


Hills 


1976) 



i s attributed to the deprival of their envelop e during close encounters with environment stars 



( jDjorgoski et al.lll99ll ; iBeer fc Daviesll2004l ). and blue-stragglers, which are main sequence 



stars more massive than the turn-off stars, may result fro m direct collisional coalescence 



and/or binary merge of two or more stars ([Leonard! Il989l . see Mapelli et al. 2006, an d 



Leigh et a l . 2007 and references the rein for recent works); see also reviews by iBailynl (119951 ). 



Hut et al.l (120031 ) and lFerrarol (120061 ). In particular, the inflation of the n umber of known blue 



stragglers boosted by the observations using Hubble Space Telescope ( Ferraro et al 



1999 



Paltrinieri et al. 2001; Ferraro et al. 2003, 2004; Sabbi et al. 2004 



Beccari et al. 



1997 



2006 



Warren et al.ll2006l ). suggests that a significant fraction of stellar populations undergo such 
encounters with neighboring stars. Recent observations with the Chandra X-ray Observatory 
indicate a link b etween the numbers of X-rays binaries and the stellar encounter rates in 
globular clusters JPoolev et al.ll2003l : IPoolev k Hutll2006h . 



Furthermore, there is a longstanding problem of large star-to-star variations in the 
surface abundances of light elements such as C, N, O, Ne, Mg and Al. Some giants in 
globular clusters exhibit the anomalous surface abundances that cannot be explain in terms 
of the nucleosynthesis and material mixi ng in the sta r s within the curr ent standard framework 
of stellar evolution (e.g., see reviews by lKraftlll994l ; iDa Costal I1997T ). Since these anomalies 
are observed only in globular clusters but not from field giants in the Galactic halo, it 
is natural to search for their origin(s) in the differences between the environment in the 
globular clusters and in the Galac tic halo, and h e nce, t o consider them as a evidence of 
the star-star interactions. In fact, iFujimoto et al.l (119991 ) have proposed a scenario for the 
formation mechanism of these abundance anomalies involving hydrogen shell flashes in red 
giants, as a result of internal mixing, triggered by the deposition of angular momentum 
into their envelopes during a close encounter with other stars. It is demonstrated that this 
extra-mixing model can reproduce the observed relationship such as correlati on and scatter 



in the anom alous abundances of Na and Al and the Mg-Al anti-correlation (lAikawa et al. 
200ll . l2004h . 



Recently, s imilar abundance ya riations are found among unevolved stars of turn-off 
and sub-giant (jGratton et al.l l200fl ). It is true that the variations and anti-correlations 
between CN and CH bands have been reported not only for giants but for stars of upper 
main sequence, which may be taken to suggest the presence of abundance anomalies i n 
unevolved stars (jSuntzeff fc Smithlll99ll ; lBriley et al.lll992uCannon et al.lll998l ; ICohenlll999l ). 
It has been argued that these facts refute the evolutionary scenario that the abundance 
anomalies are produced during the evolution along the giant branch and favor the primordial 
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scenario that the s tars were born of gas already subject to the anomalous abundances (e.g., 
Sneden et al.ll2004j ). As a possible compromise, the recycling scenarios have been proposed 
that the inhomogeneity is due to the surface pollution by accreti ng the ejecta of an omalous 
abundances from the erstwhile AGB stars of intermediate masses ( IThoul et al.ll2002l ). or that 
the second-generation stars were born from the gas pollu t ed wi th the ejecta of AGB stars 
(ID'Antonalhood ; b'Antona fc Caloil I2004T ) . Iventura et all (I200lh argue that the burning at 
the bottom of convective zone in low metallicity AG B stars can produc e the observed O- 
Na and Mg-Al anti-correlations. On the other hand, iFenner et al.l ( 2004 ) cast doubt upon 
the relevance of AGB ejecta to the observed anomalies. In addition, the scenarios involve 
serious difficulties both in the mass supply and in the overabundances of CN and s-process 
elements, attendant with the third dredge-up during the evolution; the necessary amount of 
mass ejecta only to cover and disguise the surface convection of giants may well exceed the 
total envelope mass that can be ejected from the erstwhil e AGB stars, and the enrichment 
of s-process elements is never observed (j James et al.l 12004] ). 



For the evolutionary scenario, it can also be argued that the abundance anomalies 
are printed onto the surface of unevolved stars through the mass transfer during the close 
encounters with such giants that have already developed the abundance anomalies; the 
anomalous abundances are themselves attributed to the deposition of angular momentum 
i nto the convective en velope of giants during prior close encounters with environment stars 
(IShimada et al.ll2003l ). The surface convectio n of population II ma in sequence stars contain 



mass of 3 x 10 3 M near the turn-off (e.g., see lFujimoto et al.lll995l ). and hence, the accretion 



of mass of this order may suffice to disguise their surface abundances with those of accreted 
matter. The evolutionary scenario with star-star interactions during the red giant branch 
taken into account are free from the above difficulties and have a fair prospect of giving a 
satisfactory explanation to these inhomogeneo us anomalies. Re cently, the helium production 
by this extra-mixing mechanism is discussed (ISuda et al.l 120071) in re levance to the splitting 



of main sequence branch, observed from uo Cen ( jBedin et al. 
jPiotto et al]l2007h . 



2004 ) and from NGC 2808 



One of the aims of this paper is to investigate whether the orbital angular momentum 
can be transferred into the envelope of giant from the orbital motion and whether the main 
sequence star can accrete the envelope mass from giants enough to disguise their surface 
layer with the accreted matter through star-star interactions. It is argued that the rotation- 
induced mixing requires (differen tial) rotation of ~ .01 times the local critical rate from 
the energetic viewpoint (e.g., see iFujimoto et al.lll999l ). although the proper theory is yet 
to be established. From the observations, the horizontal branch stars are known to display 
a bimodal distribution of rotation velocity with the fastest rotators at velocity f ro t sin i > 
30 km s -1 (where i is the inclination angle of spin axis) on the cooler side of horizontal 
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branc h where < 150 00° K (IPetersonl Il983l: Peterson et all 1199a ICohen & McCarthy 
19971 iBehr et all boOOal Jbl: [ Recio-Blanco et al.ll2002l ). If we neglect the angular momentum 
loss during the transition to the horizontal branch, such rapid rotations require the angular 
momentum corresponding to the rotation rate of an order of Q ~ COI^rg a t the tip of 
red giants (^rg being the critical rotation rate at the surface). This poses a problem of the 
origin of angular momentum since in the low-mass stars, the angular momentum is effectively 
extracted by magnetic br aking during the main sequence phase and by ma ss loss during the 
red giant phase (e.g., see iRecio-Blanco et al.ll2002l ; ISuda fc Fujimotoll2006l ). 



The star-star interactions have been proposed as the mechanism(s) to form the unusual 
stellar ob jects discussed abo ve and studied by many authors. Among the analytical ap- 
proaches, iFabian et al.l (119751 ) first presented an idea and evaluated the possibility that the 
low-mass X-ray binaries are produced through tidal dis sipation during the t wo bo dy encoun- 
ters involving a neutron star or low-mass black hole. iPress fc Teukolskyi (119771 ) developed 
the linear perturbation theory of the two-body tidal capture mechanism to derive a general 
formula for the amount of orbital energy, deposited into the o scillatory modes of stel lar enve- 
lope during a periastron passage. iLee fc Ostrikerl (119861 ) and McMillan et al.l (119871 ) worked 
out the cross sections for the binary formation via tidal capture of a main sequence star and 
McMillan et al.l (Il990l ) that of a red giant, respectively. These studies are, however, limited 
to the linear regime and can not deal with the non-linear effects such as the mass transfer 
between the stars and the mass loss from the stars owing to large deformations of the stars 
by tidal force. 

In order to estimate the non-linear effects during the close encounters, numerical sim- 
ulations are necessary. Among the numerical approaches, most studies have been devoted 
to understanding the resultant offspring of the stellar interactions. Simulations of tidal en- 
counters have been pe rformed for the var ious combinations of stars, e.g., a main sequence 
and a red giant star (IBenz fc Hillsl Il99ll ): a neutron star and a main sequence star or a 
red giant star in an att empt to explain th e formation of the low-mass X-ray binaries and 



the m illisecond pulsars (jDavies et al.l Il992l ; iRasio &: Shapiro I Il99ll ; iDaviesI Il995l ; ILee et al. 



19961 ): ma in sequence stars in encounter and c ollision, aiming at the formation of blue 
stragglers (ILai et al.lll993l ; lLombardi et al.ll2002l ): and a red giant star a nd a neutron star 



in re l ation to the formation of pulsars or ultra-compact X-ray binaries (IRasio & Shapiro 



199ll ; lLombardi et al.l 120061 ) . These studies have been performed exclusively by using the 
smoothed particle hydrodynamics (SPH) method, except for the encounters involving a 
massive black hole and a star, which were calculated by using a three-dimensional Euler 
hydrodynamic code fchokhlov et alill993al lbT). Recently SPH simulations are a pplied to the 



evolution of a giant planet throu gh the tidal interactions with a sun-like star (IFaber et al. 
20041 ; llvanov fc Papaloizoull2004h . 
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The former hydrodynamic simulatio ns, especially t hose o f the tidal encounters between 



a red giants and a main sequence star by lDavies et al.l (Il99ll ). have been carried out with a 



relatively small number of SPH particles and for limited range of parameters. Their results 
are thought to be subject to limitations arising from low-mass resolutions since mass involved 
in the interactions decreases as the periastron distance increases, and the SPH method may 
not give a valid description of such situations where mass scales are as small as that allotted 
to each particles. In our problems of surface pollution, we deal with the accretion of mass 
~ 10 _3 M Q . Simulations with finer mass resolutions, and hence with larger particle numbers, 
are necessary to investigate such encounters involving the transfer of mass of this order. It is 
also desirable to perform simulations for a wide range of parameters, such as the periastron 
distance, the red giant models in different evolutionary stages, and the mass of main sequence 
stars in order to obtain a realistic and general information about the characteristics of the 
tidal interactions. 

In this paper, we first carry out simulations of tidal interactions between a red giant and 
a main sequence star by using SPH method. We make a detailed analysis of the amounts of 
energy and angular momentum, transferred from the orbital motion to the oscillation and 
spin of red giant, and the amount of mass, lost from the red giant and accreted onto the 
main sequence stars. Based on the numerical experiments, we then attempt to clarify the 
parameter dependence of these characteristics and to formulate the quantitative outcome 
as simple functions of the stellar parameters and impact parameter of encounters. The 
present results are applied to investigate the relevance of the scenario that the star-star 
interactions give rise to the abundance anomalies observed among not only giants but also 
main sequence stars in globular clusters. The derived formulae will be useful to perform 
simulations of dynamical evolution of stellar systems with the effects of stellar interactions 
taken into account. 

The organization of the paper is as follows. In next section, we describe our numerical 
methods, including the set-up of the initial conditions, the models of the red giant, the 
treatment of accretion and the determination of the viscosity of red giant models. In §3, 
we present the results from our simulations with the discussion of the non-linear effects of 
tidal interactions. In §4, we derive fitting formulae for the energy and angular momentum 
deposited into the envelope of the red giants and the mass accreted onto the main sequence 
stars during the tidal encounters. The conclusions follow in §5, with the discussion about 
the application to the globular clusters. 
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METHOD OF NUMERICAL COMPUTATIONS 



In the present work, we use the three di mensional, smoothed particle hydrodynamics 
(3 D-SPH) cod e , orig inally developed by Benz (jBenzlll99Cll ; iBenz et al.lll990l ). and extended 
by iBate et al.l (119951 ) . The v ariable smoo thing length is adopted with the hierarchical tree 
method, originally written by iPresd (119861 ). to make the list of particles in the closest neigh- 
borhood of the particles. The kernel and the integrating method in our code are respectively 
the standard cubic-spline kernel and a second-order Ru nge- Kut t a- Fehlberg integrator with 
individual time steps for each particle (IBate et al.lll995l ). 



In our cases, the timescale of periastron passing is much shorter than the Helmholtz- 
Kelvin timescale in the envelope [thk — 1.5 x 10 4 (Rt {G /20R (3 )~ 1 (Lrg/ Lq) -1 y r , where .Rrg 
is the surface radius of red giant]. Accordingly we assume the adiabatic relation for the gas 
in the envelope of red giants. In actuality, the code takes into account the change of the 
entropy due to viscous dissipation, although it may have only a minor effect since we deal 
with the tidal interactions at large distance, not accompanied by large shock dissipation. 

Our SPH code uses the standard form of artificial viscosity with two free parameters 
q;sph and /3sph, which respectively control the strength of the shear and b ulk viscosity compo- 
nents and that of a second-order, von Neumann- Richtmyer-type viscosity (IMonaghan fc Gingold 
19831 ). It is known that the linear artificial vi scosity can be r e duced to the Shakura-Sunyaev 
viscosity prescription in the continuous limit; iMeglicki et al.l (119931 ) derived a relation which 
connects the viscous force with the linear artificial viscous parameter cksph- If the density 
varies on a length-scale much larger than the velocity, the shear viscosity, u, is written in 
terms of the artificial viscosity parameter asPH; i n the form: 



(l/10)a S PHcA 



where c s is the isothermal sound velocity and h is the smoothing length (jOkazaki et al.ll2002l ). 
In the envelope of red giants, on the other hand, we may relate the shear viscosity to the 
eddy- viscosity of convective motions, Veddy, evaluated at 



v, 



eddy 



•It 



(2) 



where f con v is the averaged velocity of convective elements and estimated from the mixing 
length theory with the mixing length, Z m j x . For a red giant model of mass O.8M and 
the metallicity [Fe/H] = —1.5, the eddy viscosity is found to be nearly constant around 

^eddy 

~ 1001/0 (see e.g. 



7 x 10 15 cm 2 s 



1 in the envelope when the radius ~ 20R Q and the luminosity 
Suda &: Fujimotoll2006l ). Since c s ■ h ~ 7 x 10 17 cm 2 s~ x on average for 



the red giant models constructed with SPH code (see below), we may approximate the eddy 
viscosity with a choice of «sph = 0.1. We perform the simulations with the two values of 
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linear artificial viscosity parameter, osph = 1.0 of common use and «sph = 0.1 in order to 
see the effects of viscous forces. As for the non-l inear artificial vi scosity parameter, we follow 
the usual prescription and set Psph = 2 ■ agpn ( IBate et al.lll995l ). 



2.1. Initial Conditions and Approximations 



Our simulations consist of two steps, i.e., we first make the initial models of red giants 
in hydrostatic equilibrium with SPH particles, and then, follow the encounter with a main- 
sequence star. Our red giant models are constructed with a total of 50,000 SPH particles of 
equal mass in the envelope and the core approximated by an appropriate external potential, 
while the main-sequence star is treated as a point mass. We start the encounter simulations 
by placing a red giant and a main sequence at a separation of 5-Rrg- Their relative velocity 
at this distance is calculated from the relative velocity at infinity, assumed to be Voo = 
10 km s _1 in this work, and the impact parameter. The red giants are assumed to be not 
rotating initially. We set the total mass at Mrg = 0.8 M Q and adopt two models at the 
different evolution stages, the one with the core mass M core = 0.32M Q and the surface radius 
Rrg = 2Oi? and the other with M core = O.48M and Rrq = 85i? ; the mass of one SPH 
particle is 0.96 x 10 -5 and 0.64 x 10~ 5 M M , respectiv ely. The former model is taken to have 
the same model parameters as by lBenz &: Hillsl (Il99ll ) and lDavies et al.l (Il99ll ) who use 7,132 
SPH particles of unequal masses, and the latter corresponds to the structure realized near 
the tip of the red giant branch. For the main sequence star, we take two different masses of 
M MS = 0.6 and O.8M . 



2.2. Red Giant Models 

Thee envelope structure of red giant can be reproduced by placing the envelope mass 
un der the influence of the gravi ty of core, modeled as a sphere of uniform density, according 
to iFujimoto fc Tomisakal (Il992l ). By solving the equations of hydrostatic equilibrium with 



an additional gravity g of the core, expressed as 



. -GM coie r/R 3 coie for r < R 

core > / o \ 

V ~ 1 -GM core /r 2 forr > R c ' 



u core i 



with the core radius -R CO re, we determine the density distribution of red giant envelope; see 
appendix A for detail. For the equation of state, we assume the polytrope of P = Kp( l+1 W 
with the polytropic index N = 1.5, which corresponds to the adiabatic equation of state with 
the adiabatic exponent T = 5/3; the polytropic constant K stands for the specific entropy 
of the monatomic ideal gas. 
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Figure [T] shows the density distributions in the envelope of red giants, thus obtained, 
for the models with different surface radii of Rrg — 2Oi? and R^q = 85R & . When the 
radius and the density are normalized with the surface radii Rrq and the envelope density 
p env = M env /i?^ G , two density distributions become nearly identical except inside of the 
core, which is a feature of red giant structure unless t he mass in the envelope is much 
smaller than the core mass (IFujimoto fc Tomisakal 119921 ). In this figure, we also plot the 
density distribution in t he red giant by taking the model from the evolutionary calculation 
( jSuda fc Fujimotol 120061 ) for comparison, which exemplifies that the analytic models can 
reproduce the envelope structure of red giants very well. Moreover, we compare the models 
with different core radii of -R C ore = O.O26-R and 2R Q to demonstrate that the assumed core 
radii hardly affect the structure outside the core of r > 2R Q , and in particular, in the outer 
envelope that may take a main part in tidal deformations, while the central density differs 
greatly by a factor of 4.6 x 10 5 . In our simulations, we therefore set -R cor e = 2-R to reduce 
the amount of particles injected into the innermost region. The initial models of red giants 
for the SPH simulation are constructed by distributing particles according to these envelope 
solutions, and then, by relaxing them into hydrostatic equilibrium with an artificial damping 
force on the particles. The relaxed distribution of SPH particle is also shown in the figure; 
it reproduces the structure of red-giant envelope very well except for the very surface layer 
of mass less than ~ 0.0001 M because of mass resolutions, where the variable smoothing 
length, h, ~ 0.1R RG . 



2.3. Accretion onto the Main Sequence Star 

We assume that the main-sequence star, treated as a point mass perturber, accretes any 
SPH gas particles that enter within the accretion radius, r acc , which is defined as half the 
Roche- lobe radius, Rl, calculated under the assumption of a circular orbit at a instantaneous 
separation, D, between the main sequence and the red giant, and given by; 

r acc = 0.5R L = 0.5,0(0.38 + 0.2 log q) (4) 

i n which the exp ression for Rl is valid for the mass ratio 0.3 < q = MmsA^rg < 20 



(jPaczynski I Il97ll ). A factor of 0.5 is adopted in order for the main sequence star not 
to artificially accrete unbound particles that happen to enter its Roche lobe. We have 
confirmed that the number of accreted particles is nearly the same with a smaller accretion 
radius r acc = 0.1-Rl- As for the accreted SPH particles, the mass, momentum and angular 
momentum that they carry are added to the point mass of main sequence stars. 
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3. RESULTS OF SIMULATIONS 

We have carried out 25 simulations of tidal encounter with a red giant of mass 0.8M Q 
and a main sequence star by varying the impact parameter, b, for eight sets of parameter 
combinations with the two different red giant models, the two different main sequence star 
and the two choices of artificial viscosity parameters. The model parameters are summarized 
with model identifiers in Table [TJ We adopt relatively heavy main-sequence stars of mass, 
0.8 and 0.6 M , based on the fact that the mass segregation may proceed to enhance the 
abundance of relatively massive stars in the cluster cores where the close encounters are 
expected to occur more frequently because of larger stellar density. In this table we give the 
periastron distance, r p , instead of the impact parameter b, which is given for a hyperbolic 
orbit as 

b 2 = r 2 p [l + 2G(M RG + M MS )Av4], (5) 

where Voo is the relative velocity before the encounter and set to be Voo = 10 km s" 1 in the 
present work. We also define the ratio, 77, between the critical angular velocity, I^rg, of 
rotation at the initial surface of red giant and the angular velocity, fi pass , for the circular 
orbit at the periastron distance, r p , as a measure of the closeness of encounter; 



^rg _ ( M RG \ 1/2 / r p ^ 3/2 



^pass \M RG + M MS J \Rrg 



(6) 



following iPress fc Teukolskyi (119771 ). In addition, we give the characteristic results of simu- 
lations, the energy and angular momentum transferred into the red giants from the orbital 
motions, and the masses, accreted onto the main sequence stars and lost from the systems; 
also listed are the periods, semi-major axes and eccentricities of orbital motions for the mod- 
els that yield bound systems and the models that end up with positive orbital energy are 
denoted as fly-by. 

Figures [2] gives the snapshots showing the variations of surface density, S, projected on 
the orbital plain for Model a8rgl (R RG = 20R & , M MS = O.8M , a S pH = 1.0, r p /R RG = 1.75 
or 77 = 1.64); the contours, separated by 0.2 dex, are plotted in the range of 10 — 0.001 times 
the average surface density, S env = M env /iiR RG , and open circle denotes the accretion radius 
of the main sequence star. Numerals in right-bottom corner give the elapsed time from the 
onset of simulation in units of d ynamical time- s cale, t rg = (R RG / G M cnv Y^ 2 , defined with the 



envelope mass of red giant as in iDavies et al.l (1l99ll ). On each panel open and filled squares 



mark the gas particles, initially situated on the two separate shells on the orbital plane, as 
the indicators of stellar rotation. 



As the main sequence star approaches, the tidal bulge is raised on the surface layer of 
red giant and grows toward the main sequence star. The oscillations of I = 2 f-modes are 
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predominantly excited as predicted from the linear perturbation theory. In the outer shells of 
a few 10 % in mass fraction, the deformations greatly elongated toward the perturber develop 
into non-linear regime, as seen from filled squares; in the interior, on the other hand, the 
perturbations remain small in linear regime, as observed from the location of open squares, 
and the gas almost stay at rest in the still deep interior. The outer deformations cannot keep 
pace with the motion of perturber as it is accelerated because of the timescale of passage 
of perturber comparable to that of oscillatory motions and of initially small rotation rate of 
the red giant. The lag of tide develops as the perturber approaches to the periastron. Later 
around time of 6.5 trg after the periastron passage, the gas streaming out of the red giant 
starts to accrete onto the main sequence star, now separated by D ~ 4Rrg- 

At the same time, there appears an interesting non-linear feature in the vicinity of the 
surface of the red giant star. As a result of the I = 2, /-mode oscillations, a density ditch 
is formed near the interface of upward and downward motions, as seen from the panels in 
the middle row. Figure [3] shows an enlarged picture of the velocity structure when the ditch 
is formed. Since the expanded mass elements gain larger tidal torque than the compressed 
mass elements, the former overtakes the latter while contracting to form an eddy-like velocity 
structure of counter-clockwise rotation. The rotation of Lagrangian shells and the resultant 
deposition of angular momentum into the red giant star due to the tidal torque proceed 
spectacularly in the surface region, dominated by the non-linear effects. Such features can 
no longer be the case in the linear theory, although the mass involved in the non-linear 
deformations is small, as seen from th e movem ent of open squares. The similar non-linear 



effects are reported by lKhokhlov et al.1 (Il993al ). who study tidal encounters between a poly- 



tropic star and a black hole, although they assume large periastron distance, and hence, the 
deformations remain nearly axial-symmetric without mass transfer, which is different from 
our case. 

This model illustrates an example that ends with the formation of a bound system 
after the encounter, as seen from Table [lj In Model e8rgl with the same parameter but for 
the smaller shear viscosity of q;sph = 0.1, the deformations are identical during the earlier 
phase of time ~ 5trq with those in Fig. [2], and the effect of smaller viscosity is plainly 
discernible only in deeper ditch that develops after 6 time units. Accordingly, there are only 
small differences in the results in Table [TJ This suggests that the transfer of energy and 
angular momentum is attributed solely to the phase-lag of the deformations, dynamically 
generated in the red giant envelope behind the perturber passage. The interactions to 
exchange these quantities predominantly occur near the periastron passage of the smallest 
separation while the deformations are still growing. On the other hand, the effects of viscosity 
become important only after the deformations contract to generate a strong shear near the 
stellar surface, and hence, the value of viscosity hardly affects the transfer characteristics. 
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The models of larger impact parameters such as Model c8rgl ofr] = 2.4 result in a fly-by 
encounter. The non-linear deformations in the outer shells of red giant are weaker than those 
of the model of closer encounter rj = 1.64 in Fig. [2J Because of slower angular velocity of 
perturber, the tidal bulge stretched out toward the perturber is relatively slimmer to form a 
chimney like structure. The ditch and the eddy-like structure, generated on the surface, are 
also in smaller scales. 

The encounter with a less massive perturber results in smaller transfer characteristics, 
a part of which is due to a larger periastron passing time, rj, when compared at the same 
periastron distance. Model c6rgl of Mms = 0.6M Q and r p = 2.00-Rrg (v = 2.14) exhibits 
a similar chimney like structure as Model b8rgl of the massive perturber with the same 
periastron distance, but the surface deformations are slightly smaller because of larger rj. 
The overall transfer characteristics lie between those of two massive perturber models, Model 
b8rgl with rj = 2.00 and Model c8rgl with rj = 2.40. 

The encounters with red giants at later evolutionary stages are exemplified by the mod- 
els of the larger radius, -Rrg — 85i? . When models with the same encounter closeness 
parameter rj are compared, the overall characteristics of interaction, i.e., the structure of 
mode oscillations, the development of deformations and the accreted process, are very sim- 
ilar in all models, despite of the large difference in the physical distance scales. This is 
attributable to the nature of envelope structure of red giants, i.e., to the self-similar density 
structure, as shown in Fig. [TJ The resultant transfer characteristics also turn out to be very 
similar when we subtract the effects of smaller envelope mass. 

The overall features of gas streaming are as follows. For small periastron distances as 
in these encounters, the non-linear effects in the tidal interactions are important and the 
surface density profile becomes highly asymmetric in the outer shells of red giants. As the 
perturber approaches to periastron, the tidal bulge is exited in the red giant and elongated, 
first directed to the perturber. Since it cannot catch up with the motion of the perturber 
because of initially slow rotation, gas particles from the surface of red giant chases after the 
perturber to gain the energy and angular momentum. Some of them eventually get captured 
by the gravitational potential of the perturber after periastron passage. On the other hand, 
most of gas involved in the tidal bulge falls back onto the red giant with gained angular 
momentum, which produces a non-linear feature of eddy-like structure in the surface region. 
Even in the case of binary formations, the orbit is highly eccentric, and hence, as two stars 
go away from each other the tidal bulge becomes slender. 



3.1. Non-Linear Deformations and The Evolution of Differential Rotation 



The transfer of orbital energy and angular momentum to the red giant is characterized 
by the following three time-scales; (i) the periastron passing time scale of main sequence 
star, r pass (= l/£l pass ), which is related to the variation of the external perturbing force: 
(ii) the dynamical time-scale of the envelope of red giant, r dyn (= l/f2 RG ), which is related 
to the stellar oscillations in response to the external force: (iii) the viscous time-scale in 
the envelope of red giant, r vis (= i?^ G /i/ e dd y ), which is caused by the convective eddy in 
the red giant envelope. In our case, r vis <C r pass < Td yn . For a slowly rotating red giant, 
the tidal bulge tends to fall behind the accelerated motion of perturber, and hence, a tidal 
lag is formed dynamically to carry the energy and angular momentum from the orbital 
motion into the oscillatory motions and rotation of red giant envelope. The viscosity plays 
secondary roles in the transfer of energy and angular momentum, as seen in Table [TJ The 
accreted mass depends little on the assumption of viscosity either. Furthermore, we see 
that the accreted mass takes nearly the same values regardless of the difference in the red 
giant models when compared among the models of the similar values of AE j (GM^ G / -Rrg) , 
instead of the encounter closeness parameter rj. This is indicative that the mass accretion is 
determined by the same process of energy deposition process. 

Figure H] shows the time variations of angular velocity, fi s heib relative to its Keplerian 
angular velocity, CI^q — (^MrgZ-Rrg) 1 ^ 2 ' aver aged over the gas particles in the three La- 
grangian rings on the orbital plane, initially located at the shells which contain the mass of 
red giant (including the core mass) by 95%, 90%, and 80%, respectively. The tidal torque 
excites the oscillations of f2 shell of periods ~ 2trg ; corresponding to the I — 2, m — ±2, f- 
mode. As the perturber approaches, the oscillations develop precipitously to reach the 
maximum strength near the periastron passage (at t ~ 4.5trg and 3.0trg f° r the models 
with _Rrg = 20i? Q and 85i? , respectively) with greater amplitudes and longer durations 
in outer shells. Accordingly, the turnover of rotation velocity delays in the outer shell, and 
at the same time, the mean rotation velocity of oscillation increases in prograde direction, 
reflecting the injection of angular momentum due to the tidal torque, which is stronger in 
outer shells. The amplitudes and mean- values of oscillations are greater for closer encounter. 
Between the different red giant models with the same closeness parameter 77, the model with 
the larger radius entails smaller variations in the rotation rate, normalized with respect to 
the characteristic rotation rate, fisheii/^RG! this is true even if we take into account the fact 
that the shell of the same mass fraction is 50% deeper in the envelope because of smaller 
envelope mass. It should be noted however that the net amount of angular momentum trans- 
ferred is larger, though slightly, for the red giant model of larger radius because of larger 
critical angular momentum (R^qQrq). 
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As for the effects of viscosity on the time variations of fi s hcii, the largest one appears 
in the difference in the minimum value after the periastron passage; for the smaller viscous 
parameter, it decreases to be smaller, and along with the phase delay of outer shells, the 
eddy-like structure of counter-clockwise flow becomes stronger, as stated above. The time of 
minimum f2 s heii coincides nearly with the time of the strongest eddy in Figs. [2j Stronger shear 
produced by the eddy-like structure in turn enhances the dissipation and inward transport of 
angular momentum; in panels a and a', we see the amplitudes of oscillations in the outer two 
shells get smaller in the second and later cycles for the models of smaller viscosity, though 
the overall similarity holds , in particular, in the shifts of mean-values of oscillations. 

Figure [5] illustrates the evolution of the radial distribution of angular velocity, averaged 
over the gas particles between cylinders, perpendicular to the orbital plane, with the outer 
and inner radii, separated by 0.1-Rrg- The outermost layer is first accelerated and pulled most 
outwards to run after the perturber of the angular velocity, f2 pass , at periastron passage (~ 
0.61^rg around the time ~ 6trq)- Then the deposited angular momentum is redistributed 
gradually into the interior. By ~ 40trg, the most of interior up to the radius R ~ 1.1-Rrg 
tends to rotate uniformly while there remains differential rotation in the outer expanded layer 
of lower density. For the smaller viscosity parameter of a = 0.1, the viscous process works 
slightly more slowly and a uniform rotation is reached only in the interior of R < 0.8-Rrg 
by ~ 40trg with stronger differential rotation in the outer shells. The timescale for transfer 
of angular momentum inside the red giants is shorter than that estimated from the eddy 
viscosity, the latter of which is ~ 1 yr and 10 yr for the case of ctspH — 1-0 and aspn — 0.1 
respectively, while trg — 0.02 yr. This also indicates that the non-linear effects, seen from 
Fig [3j must mainly contribute to the redistribution of the angular momentum. 

3.2. Transfer of Energy and Angular Momentum and Mass Accretion 

We evaluate the change in the orbital energy, AE, from the difference of the kinetic and 
potential energy as 



where n is the reduced mass: r and v are the relative distance and velocity, respectively, 
and the subscript denotes the quantities before the encounter. We may well assume that 
the change in the orbital energy is equal to the energy transferred to the red giant since 
the energy carried away by the particles that escape from the system is much smaller and 
escaped mass is smaller than accreted mass by nearly an order of magnitude as seen from 
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The total angular momentum, AL, transferred into the red giant is estimated by sum- 
ming up the specific angular momentum for all the particles constituting the red giant 
envelope around the core; 

AL(t) = 5^mi(r< - t-r G ) x (t% - vrg), (8) 

i 

where m{, Vi and V{ are respectively the mass of the i-th gas particle and its position and 
velocity vectors, and r RG and Wrg are the position and velocity vectors of the core of red 
giant, respectively. We exclude from the summation the gas particles which have accreted 
onto the main sequence and those which have escaped from the system; the latter particles 
are defined as satisfying the following two conditions; (1) the total energy, i.e., the sum of the 
thermal, kinetic and potential energy, of a gas particle is positive and (2) the radial velocity 
is positive when measured from the center of mass. 

Figures [6] and [7] show the time variations of orbital energy and angular momentum 
(-Eorb = f^QV^ — AE and L or b = fiobvoo — AL), respectively, for the 2Oi? models. As the 
periastron is approached, both decrease rapidly and reach the minimum after the periastron 
passage. Then, they turn to increase gradually to resume the loss up to ~ 28% at the largest 
case, and approach to asymptotic constant values. Figure [8] shows the time variation of mass, 
M acc , trapped by the perturber; the beginning of mass accretion coincides with when the 
orbital energy and angular momentum hit the minimum, and the increase in the accreted 
mass follows the curves of the latter' recovery. This indicates that the orbital energy and 
angular momentum once received by the surface matter are slowly returned back to the 
orbital motion by the accretion process. 

At the end of our simulations, the motions of two stars tend to settle in the asymptotic 
orbits, and the characteristics no longer change. We present the values of ALJ(^), AL(£ E ) 
and M acc (t^) at the end of our simulation at t = t E = 20 or 40r RG in Tabel [TJ The 
deposited energy and angular momentum into the red giant envelope increase precipitously 
with decreasing periastron distance, and for r p < 2i?R G , both of them become appreciable in 
comparison with the gravitational binding energy and the angular momentum corresponding 
to the critical rotation of red giants, respectively. The angular momentum of AL(tg) > 
0.01/rg^rg (where Jrg is the moment of inertia of red giant) can be deposited, which is 
necessary to explain such fast rotators as observed for the horizontal branch stars in the 
globular clusters. The accreted mass onto the main sequence stars also shows a similar 
tendency, amounting to be comparable with the mass in the surface convective zone of 
main sequence stars ne ar the turn-off stars in the globular clusters (~ 0.003 M ; see, e.g., 
Suda &: Fujimotoll2006l ). As the perturber mass decreases by 25% from Mms = 0.8 to O.6M , 



the deposited energy and angular momentum decrease by a factor of 1.6 — 2.2, and the accrete 
mass decreases by a slightly larger factor of 2.3 — 2.5. When the two red giant models are 
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compared at the same encounter closeness parameter 77, the interactions tend to be weaker 
for red giant models of larger radius, when normalized with respect to their radii, giving 
smaller AE(t E )/(GMl G /R RG ), AL(t E ) / (GM^R^) 1 / 2 , though the differences remain less 
than a factor of 2. This is attributable mainly to the smaller envelope mass involved in the 
tidal deformations (decreasing by 50 %) for the model of the larger radius. 

In our simulations the border between the tidal capture and the fly-by, i.e., whether 
two stars form a bound system or fly apart after the encounter, lies in the range of 2.14 < 
t p I Rrg < 2.25 for 2Oi? and 1.41 < r p / R RG < 1.65 for 85R Q , respec tively. Our tidal c aptur e 
limit is somewhat larger than obtained from the linear analysis by iMcMillan et al.l (119901 ) . 
who give the range r p / Rrg = 1.5 — 1.7 for a capture of a O.5M dwarf (see Table [2]) and 



Bailynl (119881 ). who give the range r p /R RG = 1.0 — 2.0 for a capture of a 1.4M neutron star. 



As for the dependence on the red giants models, the case of 85_R results in the fly-by for 
closer encounters than the case of 20R Q when compared between the models of the same 
periastron distance normalized by the stellar radius. This is due to the smaller binding 
energy of envelope, which directly affects t he capture condi t ion A E(t E ) > (l/2)/i f^ o , as 
already discernible in the linear analysis by McMillan et al.l (1l990l ). On the other hand, a 
larger radius causes relatively greater effects on the angular momentum deposition and mass 
accretion; for the red giants of late evolutionary stages, even the fly-by encounters can give 
the sufficient amounts of angular momentum and accreted mass to explain the fast rotators 
of HB stars and to disguise the surface of the main sequence star with accrete matter. 



4. Parameter Dependences of Transfer Characteristics and Fitting Formulae 
4.1. Tidal Energy Deposition and Angular Momentum Transfer 



The linear perturbation theory has been developed by iPress fc Teukolskyl (119771 ) and 



Lail (119971 ) to evaluate the transfer of energy and angular momentum through the dynamical 
tides; according to their results, the parameter dependences of these quantities are given 
explicitly in eqs. (IB7I) and (IB13j) for the I = 2, /-modes (see Appendix B). In order to 
separate the effects of the secondary mass, we may define AE and AL as 

AE = AE(t E )/[(GM^ G /R RG )(M MS /(M RG + M MS )) 2 } (9) 
AL = AL(t E )/[(M RG R 2 RG n RG )/(M MS /(M RG + M MS )) 2 ]. (10) 

In the linear theory, AE and AL are expressed in terms of the transfer functions, T2(f]; Q02, 
and S 2 (r); Q 02 , ^02); given in eq. (1B8[) and eq. (1B15[) in Appendix B, respectively: here u 02 
and Q02 are the frequency and overlap integral of I = 2, /-mode oscillations and for the red 
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giants, given by interpolation formulae of eqs. fIBlOl) and (1B9I) in Appnendix B as a function 
of radius. Since the I = 2 modes dominate over the oscillatory motions even in the non-linear 
regime, we may utilize these transfer functions and seek the fitting formulae of AE and AL 
as a function of rj for the given models of red giant. 

Further, the transfer characteristics obtained by the numerical simulations differ also 
with the internal structure of red giant models, as seen above. In the linear theory, the 
dependences of the internal structure are included in the transfer functions, in particular, 
through the overlap integral Qo2- in the tidal interactions via the torque, the coupling is 
given by the moment of inertia of the envelope since the core acts as an inert source of 
gravity. As a corollary, the transfer functions can be scaled with the inertia of envelope. We 
define the non-dimensional moment of inertia, J, as 

1 = J 4nr*pdr/M RG R 2 RG ~ 0.15(M cnv /M RG ), (11) 

for the red giant models. The approximation in the rightmost member in eq. (TTTT) follows 
from the similarity of the envelope structure in the red giants when normalized with respect 
to the surface radius and the envelope mass, as shown in Fig. [TJ For a main sequence model 
of polytrope N = 1.5, we have I = J 47rr 4 pdr / MusRyiB — 0>2. 

Figure [9] shows AE/(I/2) and AL/ 1 as a function of rj for the models with the different 
red giant models of radii 20 and 85i? (denoted by circles and squares, respectively) and 
with the different main sequence stars of mass 0.6 and 0.8 M (denoted by open and filled 
symbols, respectively). They form a single curve on each panel, independen t not only of 



the pe rturb er mass but also of the red giant models. The models computed by iDavies et al. 



(119911 ) and IDavies et al.1 (jl992j) are also plotted in this figure; their models of a 0.8 M Q red 
giant and 20R & give very good agreement with ours for the encounter not only with a 0.6 M & 
main sequence star (open triangles) but also with a 1.4 M & neutron star (filled triangles). 
The encounter with neutron stars result in slightly smaller energy deposition (about several 
tenths) for close encounters of rj < 2, which is attributable to larger accretion radius of 
neutron star since the accretion of larger mass onto the perturber returns a larger portion 
of energy from the outer elongated part of red giant envelope to the orbital motion. As for 
the angular momentum, the deposition is slightly smaller in their red giant models at close 
encounters of rj < 2 than in ours, which may stem from the difference in the criterion of 
mass loss particles, giving a larger mass loss to their models, or from larger radii of red giant, 
used in the normalization, with taking into account the swell of the red giants during the 
encounter. 



In the upper panel, we also plot the results for the encounter simulatio ns o f main se- 
quence stars with a neutron star and with a black hole by IDavies et al.1 (119921 ) and lLee et al. 
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( 19961) and for the en counter of a iV = 1.5 main sequence star with a black hole by 
Khokhlov et all (ll993al ). We approximate the main sequence to a polytrope of index N — 1.5 
and take I = 0.2. These models fall very closely along the same curve as our models, and 
seemingly compose a single group despite the differenc e not only i n the mass of the main se- 
quence stars but also in the mass ratio. T he models by lLee et al.l (119961 ) give slightly smaller 
values than those by iDavies et al.l (119921 ). which may stem from the differ ent criteria and 



treatm ent of the particles that accrete onto the neutron star. Th e models by lKhokhlov et al. 



(Il993al ) show slightly larger than those by IDavies et al.l (119921 ) , which may stem from the 
neglect of the accretion effect. 

This convergence may be related to the fact that the ratio, rjuj 02 , between the timescale 
of periastron passage and the timescale of envelope oscillations decreases near to the unity in 
the non-linear regime for small rj. The similar tendency that the dependence on the stellar 
models becomes weaker for smaller 77 is also discernible in eqs. (1B7|) and (1B13j) from the linear 
perturbation theory. These results for the linear regime are also plotted in this figure, and 
the comparison with the results of non-linear regime indicates that the latter effect enlarges 
the deposition of energy and angular momentum by a factor of several and up to ten, while 
it becomes saturated and slightly dwindles for smaller 77 < 2 because of the mass accretion 
onto the perturber. 

We may take advantage of the convergence in the non-linear regime to evaluate the 
transfer characteristics for other red giant models and to seek the fitting formulae that 
express AE(1;e) and AL(t^) in terms of model parameters. We define the critical rota- 
tion energy and angular momentum of model stars as _E cr it = (1/2) I{GM rg /Rrq) and 
-^crit — -^^rg-Rrg^rG' respectively, with the dependence on the moment of inertia taken 
into account, and assume the following fitting formulae that converge to the results of linear 
theory for distant encounters; 



AE(t E )/[E CTit {M P /(M RG + M P )} 2 ] 
AL(t E )/[L mt {M P /(M RG + M P )} 2 ] 



(2/I)r ] - 4 T 2 (r ] ; Q 02 , u 02 )[l + exp(a 1 r ? 2 + b lV + c{)]2) 
(l/I)r}~ 4 S 2 {r], Qo2,^o2)[l + exp(a 2 7? 2 + b 2 rj + c0]3) 



where Mp is the perturber mass. For the red giant models of various evolutionary stages, 
the transfer functions are computed with the estimates of ujq 2 and Q02 from eqs. (1B9I) and 
(IB 101) in Appendix B. For the main sequence stars, the values of ojq 2 and Q02 are taken 
from lLee fc Ostrikerl (119861 ) for a. N = 1.5 polytrope. We may determine the coefficients in 



these formulae by applying the fitting procedure with the non-linear least-square Marquardt- 
Levenberg algorithm to our data plotted in Fig. for the red giant models of radii 20 and 
85i? . The fitting curves are plotted in the figure, which converges into a unique relationship 
in the non-linear regime for the red giant models of different radii. The fitting curves for 
the models with other radii (core masses) also derived by adopting the same data in the 
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non-linear regime, as shown for the model of 4Oi? in the figure. These fitting parameters, 
thus obtained for the various red giant models, are expressed as the second order polynomials 
of radius as: 

a x = -7.8 E-05(R RG /R Q ) 2 + 9.6 E-03(R RG /R Q ) - 0.66 
h = 9.0 E-05(R RG /R Q ) 2 - 0.028(R RG /R Q ) + 3.7 

ci = -5.4 E-0A(R RG /R Q ) 2 + 0.077(R RG / R Q ) - 4.6 (14) 
a 2 = -4.6 E-06(R RG /R Q ) 2 + O.OOlO4(i2 RG /i2 ) - 0.05945 
b 2 = 2.6 E-0A(R RG /R Q ) 2 + 0M27{R RG /R Q ) + 1.02 
c 2 = 2.2 E-0A(R RG /R & ) 2 + 0.036(R RG /R Q ) - 1.42 

Furthermore, in the non-linear regime, the deposited energy and angular momentum 
during the encounters may be given simply as the functions of ry in the form: 

AE(t E ) = E crit (Mp/{M RG + M P }) 2 exp(2.718-0.761r/- 0.386r/ 2 ), (15) 
AL{t E ) = L crit (Mp/{M RG + M P }) 2 exp(3.735 - 2.23777- 0.040r/ 2 ), (16) 

which are plotted in this figure by broken lines. Note that they coincide with the expressions 
obtained above in eqs. (fl2"|) and (|T3|) for close encounters of r\ < 4 and r\ < 10, respectively. 
These formulae are applicable not only to the encounter of red giants but also to that of 
main sequence stars of arbitrary mass and radius as long as the stellar mass is small enough 
for the surface convection to develop deep enough to be approximated by a polytrope of 
N = 1.5. 



4.2. Accreted Mass onto Main Sequence Stars 

During the encounter, matter near the very surface may gain a lot of energy and angular 
momentum from the orbital motion to expand and eventually be captured by the perturber. 
The matter accreted onto the perturber returns the acquisitions back to the orbital motion, 
which may reduce the transfer of energy and angular momentum at small rj < 2, as stated 
above. As seen from Table [H the accreted mass, M acc (t£;), turns out to be nearly the same 
among the models of the similar values of AE(tE)/{GM RG / R RG ) regardless of the radius of 
red giants. This is indicative that the amount of accreted mass is related to the deposited 
energy. In this section, we study the relationship between the accreted mass and other 
physical quantities and attempt to express the accreted mass as a function of the model 
parameters. 
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In the case of a star filling the Roche lobe in a close binary, iPaczynski fc Sienkiewics 



(119721 ) argue that the mass transfer rate is related to the excess, AR, of stellar radius over the 
Roche lobe under the assumption of polytrope; the principal part of parameter dependences 
of the transfer rate (see their eq. [A21] ) is approximated by; 

GM \ NW /M RG \ N+1S fAB.\ N+ls 



M ~** A '\rT) W R \ir L ) • (17) 

where A and M are the separation and total mass of the binary system, respectively. Here 
we neglect the weak dependence on the mass ratio, and in particular, take the Roche radius 
Rl/A « 0.38 (cf. eq. [4]). Although the flow is not in steady state in our case, the timescale 
of flow through the inner Lagrangian point is slower than the dynamical timescale of stellar 
envelope, and hence, we may assume the same dependences for the accreted mass. Further 
since the orbit is eccentric and not circular, it is difficult to estimate AR/R^ exactly. And 
yet, it seems natural to AR/R^ to be related to the deposited energy and we may well 
assume the following relation; 

{d<f) L /dr)AR L AR L R RG f(AF/lrM 2 m n (m) 
GM/R RG ~ -R^l^ ~ f^E/[GM KG /R KG \) (18) 

with taking account of the work against the gravitational potential (0l) at the Roche lobe 
surface. By using the relation between the polytropic constant K and the stellar surface 
characteristics in eq. ( 1A7I) . we then have; 

M ~ M RG (^f-^ 2 f(AE/[GMl G /R RG ]f^. (19) 

iWenv * ^pass 

Here we have replaced the orbital angular velocity by the instantaneous angular velocity, 
^pass(= a/ GM/A 3 ), of circular orbit at periastron distance. Consequently, multiplying 
eq. (IT91 by the periastron passage time ~ fi~^ s and putting N = 1.5 lead us to: 

M acc ~ M RG {M enw /M RG )- l l 2 rff{AE/[GMlJR KG \) N+1 -\ (20) 

Figure [10] shows the accreted mass, M acc (t£;), divided by MRG(^env/A^RG)~ 1//2 ^ 3 , against 
AEits) I '{GM RG / R RG ). It is clearly seen that our numerical results, denoted by open and 
filled circles and squares, form a unique relationship, indifferent of the mass of perturber 
and of the red giant model; a power-low relation is discernible in the range of AE(tE) > 
3x 10~ 3 GM^ G / -Rrg- F° r smaller AE(t E ), the accreted mass tends to drop off from the power- 
law relationship because of the low mass resolution due to the limited number of gas particles; 
in actuality, only a few gas particles are accreted aroun d AE(tn) 



In this figure, we plot the results of the simulations by iDavies et all (Il99ll ). iDavies et al. 
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( 119921 ). and iLee et al.l (119961 ). and find that their values also fall very closely onto the same 
relationship for relatively large deposited energy, while the lower mass resolution in these 
earlier simulations (7500 and 9185 SPH particles, respectively) causes the deviation at larger 
AE(t E ) /GM^ G / R RG . If these results for the encounter of main sequence stars are included, 
the power relationship holds in the range of mass accretion over four orders of magnitude or 
more. 

The power-law fitting to eq. (1201) yields 

M acc /M RG = 3.5(M env /M RG )- 1 /V[A J E(t s )/(GM2 G /i? RG )] 1 - 93 , (21) 

which gives the accreted mass as a function of model parameters along with the fitting 
formula of AEitg) in eq. (|12p or in eq. (|T5|) . This relation implies that AR/Rl oc AE 1,93 ^ 3 . 
We show the comparison between this relation and the results of simulations in Figure [11] 
as a function of r] which is more useful than that of AE{t,E)- It shows a good agreement for 
the both red giant models. Since we may regard the deviations for large r?(> 2.5) as due to 
the low resolution in mass in the simulations, this gives a reasonable fitting for the accreted 
mass as a function of i] for any given set of the model parameters of encounters. 



5. CONCLUSIONS AND DISSCUSION 

We have performed the SPH simulations of tidal encounter of red giants with environ- 
ment stars and investigate the characteristics of stellar interactions for a variety of sets of 
parameters, the evolutionary stages of red giant, the mass of perturber stars and the assumed 
strength of viscosity as well as the orbital parameters of encounter. Based on our results and 
the other extent models, we discuss the dependences of interactions in the non-linear regime 
on the stellar and encounter parameters, and proposed formulae to describe the energy and 
angular momentum deposition to red giants, and the mass accretion onto the perturber stars 
in simple and convenient forms as a function of these parameters. Our main quantitative 
results £1X6 clS follows; 

1. We obtain the both amounts of energy and angular momentum, transferred from the 
orbital motion into the oscillation and rotation of red giants, during tidal encounters 
by numerical simulations. The angular momentum deposited in the red giants can be 
large enough to rotate the envelope at rate Q > 0.01f2 R c for the encounter of periastron 
distance r p /i? R c < 2.5 (or the impact parameter 6/-R R g < 15.7) and hence for such 
encounters that ended in the two stars flying apart. For still closer encounters, it 
increases to give the rotation rates significantly exceeding Q ~ 0.1f2 R Q. Accordingly, 
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the tidal encounter works as the source of angular momentum necessary to trigger 
rotational mixing in the red giants and also to explain the origin of fast rotators 
observed among the horizontal branch stars. For larger radius of red giant model, 
and hence, later stage of evolution, the transferred angular momentum increases while 
the energy deposition decreases since the transferred quantities are scaled with the 
stellar parameters. The fitting formulae are derived to describe these quantities as a 
function of the mass and radius of red giants, subject to the perturbation, the mass 
of main sequence stars as the perturber, and the impact parameter. Further we show 
that these transferred quantities, when normalized with respect to the momentum of 
inertia of models, are given solely as functions of the encounter closeness parameter 
77, Our fitting formulae agree well with the results of encounter simulations by other 
authors, and can even reproduce the results for the encounters of main sequence stars 
approximated to a polytrope of index N = 1.5. 

2. With aid of fine mass resolution, we demonstrate that the main sequence stars can 
capture gas from the red giant envelope sufficiently to disguise their surface with ac- 
creted matter even for the encounters that ended in fly- by with an evolved red giant. 
The accreted mass onto the main sequence as a perturber during the tidal encounters 
is shown to be in a direct relationship with the energy deposition into the red giants. 
We also derive the formula, which predicts the accreted mass as a function of impact 
parameters for given stellar parameters and are applicable to the encounters involving 
not only the red giants but also the main sequence stars. 

The derived formulae are useful in determining the periastron distance of the tidal 
capture limit for the encounter of various model parameters. They also be useful in inquiring 
whether some stellar objects in the globular clusters, for example the red giants and the main 
sequence stars with abundance anomalies and the fast rotating horizontal brunch stars, which 
cannot be explained through the framework of the normal stellar evolution, can be produced 
through the stellar interactions. In our computations, the mass accretion rate may exceed 
the Eddington limit (Mgdd — 47rci?Ms/ K e) on the surface of main sequence star for very close 
encounters, but since it remains below the Eddington limit at the accretion radius, we may 
well assume that the accreted mass mostly settles on the surface of main sequence stars with 
loosing thermal energy. In the following, we discuss the application of these formulae and 
the possibility that such stellar objects have their origins in the stellar encounters. 
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5.1. Tidal Capture Limit and Comparisons with Other Works 



In our simulations, the tidal capture limits, r p cap in periastron distance and r] cajp , in 77, 
are estimated from the condition that AEfa) = (l/2)/xVoo(10 km s) 2 with use of eq. (}T2"]) or 
eq. f|T5|) . Our estimates are given in Table [2] for the two red giant models with the perturber 
masses of 0.6M Q , O.8M and 1.4M . For the red giant of larger radius, the tidal capture 
limits decrease slightly (~ 10%) when normalized with respect to the radius of red giant, but 
increase in the physical dimensions nearly in proportion to the surface radius; they slightly 
increase with the mass of perturbe r. For the encounte r of M rq = O.8M red giant with a 
Mms — O.6M main sequence star, iDavies et ajj (119911 . Il992l ) give the periastron distances 
for the tidal capture limit in the range of 2.00 < r PiCa p/_R RG < 2.25, and our estimate, 
f>,cap/-RRG = 2.1, as seen from Table [2], resides in their range. On the other hand, our 
estimates turn out to be larger by ~ 20% than those obtained from the linear analysis by 
McMillan et all Jl99oh . as listed in the table. 



On the other hand, iKhokhlov et al.l (Il993al ) and iKhokhlov et al.l ( 11993bl ) compute the 
encounter of a polytrope star of mass M* = O.8M with a black hole of mass Mb 3> M* for 
various values of polytropic index (the relative velocity is 100 km/s at infinity). The tidal 
capture limit decreases from r] cap = 2.75 for N = 1.5 to 2.20 and 1.55 for iV = 2.0 and 3.0, 
respectively. This demonstrates that for larger polytrope index, the energy deposition rate 
of the star becomes smaller, because of the increase in the mass concentration toward the 
center and of decrease in the moment of inertia for a given mass and radius. 



5.2. Relevance to the Origin of Stars with Anomalous Abundances in 

Globular Clusters 

We may apply our fitting formulae to examine the possibility that the abundance anoma- 
lies observed both for red giants and main sequence stars in some globular clusters can be 
explained in terms of the stellar interactions during the close encounters. As the origin 
of these objects, we propose the following scenario; 1) The anomalies of red giants are 
generated through the flash-assisted deep mixing mechanism that is triggered by the injec- 
tion of angular mome ntum into their envelope during the encounter with environment stars 



(IFujimoto et al.lll999l ): 2) The main sequence stars gain the abundance anomalies as a result 
of the surface pollution by accreting matter from the red giants which have already devel- 
oped these anomalies. We evaluate whether these work under the conditions prevailing in 
globular clusters. 



First as for the point 1), iFujimoto et al 



(119991 ) argue that the angular momentum of 
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ALRc/^crit > 1/100 is necessary to induce the flash-assisted deep mixing, and the trans- 
fer of angular momentum of this order occurs during the encounter of i] < 3, and hence, 
r p ~ 2.6-Rrq from our formulae. The transfer of angular momentum may have relevance 
to the a bimodal distribution of rotation velocity that horizontal branch stars display with 
the fastest rotators distributed on the cooler (redder) sid e of the branch whereas with the 
slower rotators spread over wider range on the branch. ISuda fc Fujimotol (120061 ) suggest 
that the different modes of helium mixing mechanism may result according to when the 
stars undergo the close encounter and the deposition of angular momentum on the RGB 
and influence the horizontal branch morphology; the injection of angular momentum may 
invoke hydrodynamical instabilities due to differential rotation and invokes turbulent mixing 
to trigger the hydrogen-flash driven deep mixing, and the resultant helium enrichment in the 
envelope accelerates the evolution of RGB. The stars that experience the helium enrichment 
at an earlier epoch on the RGB have a smaller mass of helium core and hence, located on 
redder side of HB and those with later mixing epoch shifts to blue-ward: on the other hand, 
if the stars experience the close encounter near the tip of RGB, the helium-flash driven mix- 
ing, rather than the hydrogen-flash driven deep mixing, takes places and causes the largest 
decrease of helium core so that the stars are situated at the reddest end on the HB. If the 
stars experience close encounter at an earlier stage of RGB, then, they become slow rotators, 
due to angular momentum loss through mass loss on the RGB, and hence, settle on redder 
side of horizontal branch, and the stars, if experiencing it at later stage of RGB, become 
faster rotators and settle on bluer side. Finally, the stars, which undergo close encounters 
very close to the tip of RGB, become the fastest HB rotators, located on redder-most side 
of the branch. The close encounters at rj > 3 can explain the fastest rotation rates observed 
from HB stars of Q ~ O.I^k if the angular momentum is conserved during the contraction 
from the RGB to the HB. While the pristine angular momentum is lost effectively for such 
low mass stars, such fast rotation as observed for HB stars may be expected also from the 
synchronization of red giants in the binary systems of separations (more than a several au), 
and yet, it is difficult for such binaries to survive without suffering encounters in the dense 
stellar environment of globular clusters (see below eq. [22]). In other word, we may take the 
existence of HB stars of these fastest rotations as an evidence that such close encounters as 
rj < 3 take place in these clusters. 

Next as for the point 2), since the mass in the surface convective zone is of ~ 3 x 1O _3 M , 
the accreted mass of an order of ~ 1O~ 3 M suffices to disguise the surface abundances 
with those transferred from the red giant envelope with anomalous abundances. From the 
present results, it is possible to estimate the range of periastron distance, rj, that can allow 
the accreted mass of this order at r\ — 1.8 ~ 2.2 and 1.5 ~ 2.0, which correspond to 
r v = 1.9 ~ 2.1i?RG and 1.7 ~ 2.0-Rrg, for -Rrg = 2O_R and 85_R , respectively, for masses 



-24- 



range of main sequence stars of 0.6M Q and 0.8M Q . 

The timescale of tidal interactions in the environment where the stellar density is rif pc~ 3 
and the velocity dispersion v go, is estimated (with the gravitational focusing taken into 
account, because of low velocity dispersion of environment stars in the core) at: 

Tm „ 7 . 4 x ^^(^(^K^i^-y (22) 

Uf r p 10km s 2M Q 

On this basis, the timescales of tidal encounters, which can bring about the mass accretion of 
a order of 1O~ 3 M and the angular momentum transfer of ALrg / Lent <i 1/100, are estimated 
at ~ 1.5 • 10 10 yr for 20R & and ~ 3 • 10 9 yr for 85i? in the environment of n/ = 10 4 pc~ 3 
and Voo = 10 km s _1 . These time-scales seem to be too large to explain the observed these 
objects by tidal interactions as compared with the lifetimes on the corresponding stages of 
red giants, ~ 10 s yr for 2Oi? and ~ 10 7 yr for 85i? , and the accumulated number of 
close encounters, attendant with the abundance anomalies, seems to be no more than a few, 
too small to explain the observations even in rough estimates. In order for the encounters 
to be viable, these time-scales have to be shorter more than an order of magnitude, and 
accordingly, there need some mechanism(s) to enhance the frequency of tidal encounter in 
globular clusters. 



Sugimotol (119961 ) points out the importance of mass segregation in modeling dynamical 
evolution of star clusters to explain the observable number of mill-second pulsars in 47 Tuc; 
without the mass segregation, theoretical estimation indicates that the formation probability 
of binary with a neutron star is higher in u Cen having no collapsed core than 47 Tuc having 
collapsed core, although in actuality, the former has not been reported to contain any pulsars. 
Moreover, mass segregation pro motes the core collapse and m ake it more rapidly than in the 



case of single-mass component. iPortegies Zwart et al.l (120011 ) confirm that mass segregation 



enriches the core of star cluster in giants and white dwarfs by N-body simulations for an open 
star cluster. In order to explain the observable number of giants with abundance anomalies, 
about a half of total number of giants have to experience tidal encounters with field stars. It 
is necessary to see whether the mass segregation can gather most of giants, though not all, in 
the core and can increase the two-body encounter rate by more than an order of magnitude. 
In addition to the ma ss segregation, the grav o -therm al oscillations of star clusters, which 



has been prop osed by iBettwieser fc Sugimotol (119841 ) and confirmed by N-body simulation 



Makinol (119961 ). may influence the rate; since the interactions are expected to occur in the 
evolutional stage near the high density peak, it is necessary to investigate in a correct manner 
how deeply and how long such a high density state reaches and lasts to affect the encounter 
rates. 



For the encounter with the red giants of later stages, the pollution of main sequence 
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stars are possible even when the two stars fly apart after the encounters. For the encounter 
that ends in the formation of binary, it depends on the fate of two stars whether the abun- 
dance anomalies imprinted onto the companions are observable or not, and it is necessary 
to pursue the details of evolution of tidally captured binaries. There are many effects to 
influence the binary evolution, such as the spin-up and mode damping rates of red giant 
afte r the encounter , the mass transfer to the companion at subsequent periastron passages 



ex. 



Sepinsky et al.l 120071 ). the mass loss from the bound system, and also the expansion of 
red giant as it ascends the red giant branch. In addition, we should also take into account 
the tidal interactions of formed binaries with the environment stars. It is conceivable as the 
destination of tidally captured binary that either one component is liberated through the ex- 
change encounter with a third body, or the two stars eventually coalesce as a result of Roche 
lobe overflow. The possibility of the exchange event depends on the evolution timescale of 
giants and the encounter timescale; since the exchange encounter rate is propo rtional nearly 



toa semi- major axis of a binary under the assumption of point mass limit (IHeggie et al. 



19961 ). then the former is expected to occur more frequently. Accordingly, if tidal capture 
of two body encounters can contribute the modification of stellar populations, then, the ex- 
change encounters follow at larger rates and make greater contributions. This increase in the 
encounter rates with red giant in the binary may affect the statistics of the main sequence 
stars with the abundance anomalies. Further, through the exchange encounters, the red 
giants and horizontal branch stars, now losing the mass due to the mass loss and becoming 
lighter than the heaviest main sequence stars, can be ejected from the binary systems to fly 
apart as single stars. 

In summary, the enhancement of tidal interactions, necessary to explain the observed 
abundance anomalies, is expected to be provided by the formation of high density core due 
to the gravo-thermal oscillations and by the mass segregations which enlarge the fractions 
of stars of the upper mass end in the core. The proper understandings of these effects wait 
for N-body simulations of globular clusters with the stars of multi-mass spectra taken into 
account since star-star interactions as studied in the present work will play a critical role. 
In these studies of dynamical evolution of stellar systems, the formulae for the transfer of 
energy, angular momentum, and accreted mass derived in the present work serve the purpose 
for incorporation these effects into the simulations. It is also necessary to pursue the binary 
evolution and accurately explore the fate of red giants and subsequent horizontal branch 
stars with the interactions with the environment stars taken into account. 



This paper is based on one of the author' s(S. Y.) dissertation submitted to Hokkaido 
University, in partial fulfillment of the requirement for the doctorate. This work has been 
partially supported by Grant-in-Aid for Scientific Research (15204010, 16540213, 18104003), 
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from Japan Society of the Promotion of Science. 



A. The Polytrope Model of Red Giant Structure 



Fujimoto fc Tomisakal (119921 ) show that the structure of red giant can be modeled by 



a combination of two polytropes with the cool and hot components corresponding to the 
core and envelope, respectively; In particular, we may replace the cool component as a 
sphere of uniform density (i.e., a polytrope of index A" = oo) since we are interested only 
in the envelope structure, and in this case, the structure of hot components ensues from the 
following equations; 

—L = 4vrr 2 p, (Al) 
dr 

1 dP GM r 

-~T = -+9, A2 

p dr r z 

with the contribution of core gravity, g, in eq. ([3]) taken into account. With a given envelope 
mass, M env , and the surface radius, R, we may introduce the dimensionless variables as 

M r -\h,nr- r = R£, p = p 6 N , P = P 6 N+ \ (A3) 

where po and Pr, are the density and pressure coefficients, related to the envelope mass and 
radius as; 

Po = M em /AnR\ P = GMiJ4ir(N + 1)R\ (A4) 

By using these variables, we may rewrite the equations in the non-dimensional form, corre- 
sponding to the Lane-Emden equation, as 

= -0" farf^e™ (A5) 

where ^ core = R core /R and y^core = M core /M env . In the above equations, we take into account 
the hot component in the core, but its contribution to the envelope mass is negligible because 
of small core radius. 

We may obtain the envelope structure of red giants by solving eqs. (1A5I) for a given set 
of core radius and mass, (£ C ore, V^core) with the boundary conditions, 

dO/d^ = 0, <p = at f = : and B = 0, (p = 1 at £ = I. (A6) 
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The several solutions are shown in Fig. [XJ, and we see that the structure, exterior to the core, 
resemble each other as long as the core radius is sufficiently smaller than the stellar radius 
(£core 1) and unless the envelope mass is much smaller or much larger than the core mass. 
In particular, we may have relation between the polytropic constant (or entropy) and the 
surface characteristics of red giant analogous to the single polytrope; 

K = Po/pI +1/N = &^M£'»&»-\ (A7) 



B. Model Dependences of Transfer of Energy and Angular Momentum due to 

Linear Dynamical Tide 

The linear perturbation theory of dynamical tides has been developed to estimate the 
deposi tion of energy and a ngula r momentum from the orbital motion to the stellar oscilla- 
tions. iPress fc Teukolskvl (119771 ) derive an analytic formula for the energy deposition, and 



later, lLail (119971 ) extends it with the stellar rotation taken into account to give the general 
formulae for both the energy deposition and angular momentum transfer, which are appli- 
cable for rj 3> 1. We are here concerned with the red giants in negligible rotation initially. 
Furthermore, we ma y well retain only a le ading term of I = 2, /-modes, which dominate the 
dynamical tides (see lLee fc Ostrikerlll986l ). 



Press fc Teukolskvl (119771 ) give the energy loss, AE, of the orbital motion into the adi- 



abatic, non-radial oscillations for a star of mass, Mi, and radius, Ri, during an encounter of 
periastron distance r p with a point object of mass M 2 , in the following formua; 



AE 



GMj M 2 2 
Ri 



1=2,3.. 



,R\ . 
TV' 



21+2 



Ti(v), 



where the dimensionless transfer function, XJ, is defined by 



(Bl) 



Ti(v) = 2vr 2 ^ | Q n i | 2 l*> 



nlm 



m=— I 



with the overlap integral, Q n \, given by 



Q nl = I r^Cr/iZO'-^rniCOrJ, 

'0 



and with the integral, K n]xa , along the trajectory, given by 

A'rilm 



!rim 2 3 /Vi m (w 



2tt 



nil 



:b2) 



(B3) 



(B4) 
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4n(y) = / dx(l + x 2 y l cos[2 1/2 y(x + x 3 /3) + 2m tan -1 a;], (B5) 
Jo 

W lm = (-)^)/ 2 [^(/-m)!(/ + m)!] 1 /V[2'(^)!(i±^)!]. (B6) 
Here £ R D i and ^ s nl are the radial and poloidal normal mode components of the Lagrangian 

— 1/2 

displacements from the unperturbed spherically symmetric state in units of M 1 , respec- 
tively, and the symbol (— ) k in (1B6|) is to be interpreted as (— l) fc when k is an integer, while 
zero when k is not an integer. Since Ti{rj) includes the stellar mode frequency, u n i, and the 
overlap integrals, Q n \, both determined from the normal mode structures, and hence, AE is 
dependent on the mode oscillation structure of the star. 

For the leading term of I = 2, /-mode, we have the following form by separating the 
dependences on the mass and radius of stars as, 

AE 2 = {GM 2 1 /R 1 ){M 2 /M 1 + M 2 ) 2 V ~ 4 T 2 ( V : Qo2,^ 02 ), (B7) 

Air 3 
T 2 (r] : Qo2,^o2) = — | Q02 | 2 [/2o(^02) 2 + -{^Om^) 2 + h-ziv^) 2 }]- (B8) 

For the red-giants of M 1 = M RG = 0.8 M Q , the fr equency, o^?, an d the overlap integral, 



Qo 2 , of I = 2, /-mode oscillations are obtained by McMillan et al.l (1l990l ). who solve the 
perturbation equations to find the overlap integral and stellar mode frequency. We may 
evaluate the values of u 02 and Q 02 by interpolating their results as a function of the radius 
in the following forms; 

lu 02 = -1.56E-O4( J R RG / J R ) 2 + O.O112(iWit! ) + 1-84. (B9) 
Q02 = 2.29E-05( J R RG /i? o ) 2 - 0.00352(^/^0) + 0.331. (BIO) 



As for the angular momentum, AL, transferred from the orbit to the spin of primary 
star, Lai(1997) derives a general form with the effects of stellar rotation taken into account. 
In the limit of negligible initial rotation rate, it reduces to 




(Bll) 



where the transfer function, Si, is defined as 

1 

S^) = -27T 2 I Qnl I 2 "nt Yl m I K ^ I' ( B12 ) 

n m=—l 

with the overlap integral, Q n \, and the trajectory integral, K n \ m , defined above. 
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Then for the 1 = 2, f-mode, we have 

AL 2 = \jG M tR^ M M + M ?r A S 2 ( V : Qo2,wte), (B13) 

where 

(B14) 

12n 

S^V ■ Qo2,^02) = -(-=-) I Q02 | 2 ^ 2 1 { / 22(^ 02 ) 2 - h-2 (W) 2 ) 2 }] • (B15) 
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Table 1. Parameters and Characteristic Results of Encounter Simulations. 
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Table 2: Tidal Capture Limits: 

First and second columns are main sequence mass and red ginat radius, third and forth ones 
are 77 and r p at tidal capture limit for parameters of first and second columns estimated from 
our simulations and last two columns are same as that of third and forth ones but estimated 
from linear theory. 
M MS (M ) R rg (R q ) r? cap r Pc jR RG ) 3|F r^(R RG ) 
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Fig. 1. — The distribution of density in the envelope of red giants, normalized with respect 
to the envelope mass and radius (M env = R^q = 1) as a function of normalized radius. Solid 
and long broken curves denote the models of the same core radius of -R C ore = 2_R with the 
same total mass O.8M , the different surface radius R^g = 2Oi? and 85R Q and the different 
core mass M core = O.32M and O.48M , respectively, while broken curve denotes the model 
of a same surface radius R^g = 2Oi? and a same core mass M core = O.32M with the same 
total mass O.8M and the different core radius -R cor e = O.O26i? , for comparison. Dash- 
dotted curve denotes a stellar model of i?RG = 19.8i? and M core = O.362M , taken from 
th e evolutionary calculatio n of a star with mass 0.8 M & and the metallicity [Fe/H] = —1.5 
by ISuda fc Fujimotd (120061 ). Asterisks represent the density distribution of the SPH model 
with Rjig — 2Oi? for the radii at intervals of 2R Q with two additional ones near the center 
and surface. 
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Fig. 2.— Variations in the surface density, projected on the orbital plane, during the 
encounter for case a8rgl (Rrq = 20R Q , r\ = 1.64, M MS = O.8M , a = 1.0). Each panel 
shows a snapshot of contour lines at intervals of 0.2 dex in the logarithmic scale over the 
range of 10-0.001 times the average surface density, S cnv = M em /tiR^q. The origin of 
coordinate is set at the center of mass and numerals in the right bottom corner indicate the 
time in units of dynamical timescale, r RG = \J R\ G / GM erw . Filled and open squares indicate 
representative SPH particles initially located on the shells that contain 95% and 70% of total 
mass including the core mass in the interior, respectively. 
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Fig. 3. — The velocity vectors, plotted on the density contour map of panel at time 7.0 r RG 
in Fig. H 
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Fig. 4. — Time variations of angular velocity, fi s hcii, averaged over particles in Lagrangian 
rings on z = orbital plain, initially located at the shells, the interior of which contains 
95%, 90% and 80% mass of red giant (including the core mass); the vertical axis is the 
angular velocity normalized with respect to the Keplerian value at the initial surface (Qrg — 
a/ GMrg / -Rrg ) an d the horizontal axis is time in units of trg = 1/^rg)- Leftmost panel 
(a) for Model a8rgl: the second and the third panels (a' and b) for Models e8rgl and b8rgl, 
which differ from the first one in the viscous parameter and in the encounter closeness 
parameter, respectively; rightmost panel (c) for Model b8rg2 with the red giant model of 
larger surface radius. 
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Fig. 5. — The evolution of radial distribution of angular velocity, averaged over the particles 
that reside temporally between the cylinders, separated by 0.1-Rrg on the orbital plane, for 
model b8rgl; f2 RG (= a/GM rg /_R^ g ) is the critical rotation velocity at the initial surface 
of red giant, and the time is designated in the box in units of dynamical time of r RG = 
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Fig. 6. — Time variations in the orbital energy, E or ^ = (l/2)/i f^ o — AE; filled and open 
symbols denote the models of q;sph = 1-0 and q;sph = 0.1, respectively. 
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Fig. 7. — Time variations in the orbital angular momentum, (L orb = /lobv^ ~ AL RG )- Filled 
and open symbols denote the models of «sph = 1-0 and «sph = 0.1, respectively. 
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Fig. 8. — Time variations of the accreted mass, M acc , onto the secondary point mass; Symbols 
have the same meanings as in Fig. [7J 
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Fig. 9. — The deposited energy and angular momentum, normalized by the inertia, AE/^I 
(top panel) and AL/I (bottom panel) from the orbital to stellar internal motions are plotted 
as a function of the encounter closeness parameter r\ for the red giant model of radius, 20 
and 85i? (circles and squares) with the main sequence perturber of mass, 0.6 and 0.8 M 
(open and filled symbols). Also plotted are the results by other authors: a red giant of 
mass 0.8 M G and radius 20R & wit h a main seq u ence of mas s, 0.4 M & and 0.6 M , and with 



a neutron star of mass 1.4 M & by iDavies et ajj ( 119911 . Il992l ) (inverted filled, o pen and filled 



triangles): a 0.8 M main sequence with a neutron star of mass 1.4 M , by IDavies et al. 



(119921 ) (inverted open triangles) : a 0.2 M , 0.3 M , 0.5 M and 0.7 M main sequence star 



1 





1e-07 

0.0001 0.001 



AE(tE)/(GM| R /R RB ) 



Fig. 10. — The accreted mass, multiplied by the root of the mass fraction of envelope of red 
giant and divided by third power of encounter closeness parameter, rj, as a function of the 
deposited energy, AE(1;e), into the red giant envelope. Symbols are the same as Figure [9] 
and solid line represents a power-law fitting. See text for details. 
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Fig. 11.— The accreted mass, multiplied by {M env /M 1 ) 1 / 2 [2/I[(M 1 M 2 )/M 2 \ 2 } im , as a func- 
tion of 77. Symbols are same as Figure [TUJ Solid and dashed curve are derived from eq. (T2T1) 
and the fitting formula eq. (|T2l) for the red giants models of 20R Q and 85_R , respectively 
(see text for details). 



